Abstract. We find that the ROSAT source 
1RXSJ105010. 3-140431 is a cataclysmic variable with 
orbital period of 88.6 minutes and a spectrum closely 
resembling WZ Sge. In particular, emission lines are 
flanked by Stark-broadened absorption wings probably 
originating in the photosphere of a compact object. The 
Balmer absorption lines can be modeled by the spectrum 
of a DA white dwarf with 13000 < T^ff < 24000 K. The 
strong absorption lines allowed us to obtain direct radial 
velocities of the white dwarf using the cross-correlation 
technique. We find an extremely low white dwarf radial 
^^ ' velocity half amplitude, K^d = 4 ± 1 km s^^. This is 
^^ , consistent with the upper limit obtained from the Ha 
^^ ■ emission line wing K < 20 km s~^. The corresponding 
mass function is incompatible with a main sequence 
secondary, but is compatible with a post orbital period 
minimum cataclysmic variable with a brown dwarf-like 
secondary. The formal solution gives a secondary mass 
of 10-20 jovian masses. Doppler maps for the emission 
lines and the hypothesis of black-body emission indicate 
J^ , a steady state (T ~ r"'^/*) accretion disk mainly emitting 
^— V ' in Ha and an optically thicker hotspot with a strong 
p\J . contribution to the higher order Balmer lines and He I 
t^^ ' 5875. As in other long cycle length dwarf novae, evidence 
^^ , for inner disk removal is found from the analysis of the 
emission lines. 
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1. Introduction 

IRXS J105010.3-140431 (02000= +10 50 10.4, (52ooo = -14 
04 32.0) is a poorly-studied ROSAT source detected with 
a count rate of 0.06018 ± 0.01495 counts/sec (Voges et 
al. 1996). In this paper we present the first spectroscopic 
and photometric study of this object, revealing a probable 
classification as a dwarf nova in the very late stages of the 
evolutionary track of cataclysmic variables (CVs). 

2. Observations and Data Reduction 

2.1. Photometric observations 

We observed IRXS J105010.3-140431 during three nights 
on February 17, 19 and 20, 2000 with the Danish-1.54m 
telescope of the ESO La Silla Observatory. The Dan- 
ish Faint Object Spectrograph and Camera DFOSC was 
used with the backside illuminated Loral/Lesser chip. This 
CCD has 2052 x 2052 15/im pixels spanning a field of 
view of 13.'3 x 13f3. A windowing of 5 'x 5 'was selected 
around the center to ensure fast readout of the CCD. We 
imaged the field of the variable continuously for approxi- 
mately ~ 2 hours every night with the V filter and then 
measured the instrumental magnitudes of the variable and 
comparison stars using the aperture photometry facility 
in the software MIRA APR. Differential photometry was 
performed using nearby isolated stars in the field. This 
technique is robust against small changes in atmospheric 
transparency and today it is routinely used for monitor- 
ing variable stars (i.e. Howell et al. 1988). A finding chart 
indicating the comparison stars is shown in Fig. 1. The C2 
—CI light curve had a mean value of -1.04 mag and rms 
0.017 mag whereas the C3 — C2 light curve showed an av- 
erage of 1.45 mag and rms 0.020 mag. We did not use 
the bright star close to the variable as a comparison since 
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it was saturated in our frames. As C2 has a magnitude 
similar to the variable, we deduce an upper limit error of 
0.017 mag for the V — C2 curve. This light curve, along 
with the errors, arc shown in Fig. 2, suggesting that the 
star is marginally variable every night at the 2 a level, but 
the mean magnitude changes from night to night by 0.05 
mag. 

2.2. Spectroscopic observations 

Spectra were obtained on the 2.2 meter telescope at ESO's 
La Silla Observatory, on January 4-7, 1997. The EF0SC2 
and CCD # 24, along with grism # 4 and a slit width of 
1 arcsecond, yielded a spectral resolution of 3.4 A, as de- 
termined by the slit width, plate scale and pixel size, and 
a wavelength range of 4020-7470 A. He-Ar comparison 
spectra were taken every five 480-s length science expo- 
sures. Bias and flat images were taken daily. Observations 
are summarized in Table 1. 

Data reduction was performed in the usual manner 
using the IRAF software packageg. All of the two dimen- 
sional raw frames were de-biased and flat-fielded using the 
average bias and normalized flat images, and then one di- 
mensional spectra were extracted. About 30 He-Ar emis- 
sion lines provided spectral calibration functions with typ- 
ical rms of 0.4 A (18 km s^^ at Ha). Flux calibration was 
performed with the standard stars LTT 1020 and EG 21 
( |Hamuy et al. 1992| , pamuy et al. 19941 ). 



3. The analysis of the emission line spectrum 

The top of Fig. 3 shows the averaged spectrum of the ob- 
ject, characterized by Hi lines and Hei emission. Flanking 
the strong Balmer emission lines we observe the pressure 
broadened absorption profiles typical of a white dwarf. We 
are confident that these lines come from a white dwarf and 
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Table 1. The spectroscopic observing log indicating start- 
ing and ending times for the time series and nightly mean 
V magnitude. The number of science frames per night is 
also given. 



Date(UT) N HJDstart HJD^ 



V 



05/01/97 


13 


53.7688 


53.8676 


17.63(09) 


06/01/97 


14 


54.7928 


54.8188 


17.82(10) 


07/01/97 


14 


55.7700 


55.8647 


17.65(05) 


08/01/97 


17 


56.7590 


56.8704 


17.48(16) 



Table 2. Mean spectroscopic quantities. Positive equiva- 
lent widths correspond to the absorption component, ex- 
cluding the emission core. 



Line flux (erg cm 



-' S-' A- 



^) Wa (A) 



R5 


5.9E-16 


-8 


H7 


6.0E-16 


-13/+33 


He I 4470 


5.5E-16 


-1 


H/3 


6.5E-16 


-33/+38 


He I 5875 


3.4E-16 


-10 


Ha 


7.9E-16 


-129 


He I 6678 


2.0E-16 


-2 



not from a hot disk, since the helium lines are not observed 
in absorption, like in outbursting dwarf novae (for a review 
of spectra during outbursts see Warner 1995, see also sec- 
tion 4). The Ha emission is double peaked (see bottom 
panel of Fig. 3) with an equivalent width of — llOA and 
half peak separation 560 km s~^. The full width at half 
maximum FWHM in Ha is 1811 km s~^ and the full 
width at zero intensity FWZI is 2950 km s^^. These 
values are typical for dwarf nova accretion disks at qui- 
escence. Spectroscopic quantities for the averaged spec- 
trum are given in Table 2. The steep emission decrement 
suggests an origin in an optically thin accretion disk. Al- 
though no outburst has been observed yet for this object, 
it is likely that it corresponds to a dwarf nova of the SU 
UMa subtype. These systems, reviewed by Warner (1995, 
2000), consist of a white dwarf accreting material from 
a cool red dwarf at very low transfer rates, producing 
thermal instabilities (outbursts) in an accretion disk with 
recurrence times from weeks to several years. The faint 
continuum produced in the disk is hidden in the brighter 
white dwarf spectrum even in optical wavelengths, produc- 
ing the absorption wings seen around the Balmer emission 
cores. In this regard the spectrum is quite similar to that 
of the dwarf nova GW Lib (Szkody et al. 2000). 

3.1. The orbital period and ephemeris 

Derivation of stellar parameters for the binary components 
in a dwarf nova requires a precise determination of the 
radial- velocity curve of at least one of the stars. This task 
is not always possible, because in most systems the white 
dwarf and secondary star spectrum is hidden in the optical 



spectrum of the bright accretion disk. In fact, the emis- 
sion lines originating in the disk usually do not reflect the 
white dwarf binary motion because they are subject to in- 
trinsic velocity fields besides the gas rotation in the disk 
itself (e.g. Robinson 1992). This yields an intrinsic uncer- 
tainty in the spectroscopic ephemeris obtained for most 
dwarf novae. The problem has been overcome in long or- 
bital period dwarf novae by observing in the infrared the 
lines of the secondary star. In these systems the relative 
flux contribution and features of the secondary star allow 
a straightforward measurement (e.g. Friend et al. 1990). 
However, short orbital period dwarf novae have secon- 
daries too faint to be detected even in the IR with current 
instrumentation. In addition, it is difficult to derive pre- 
cise radial velocities (R.V.) using molecular bandprofiles. 
Fortunately, the disk in IRXS J105010.3-140431 is faint 
enough to reveal the white dwarf photospheric lines, so 
we have a unique opportunity to obtain white dwarf ra- 
dial velocities and the unbiased ephemeris for the system. 
The spectra were normalized by dividing them with a low- 
order polynomial continuum fit before we performed the 
radial velocity analysis. 

We traced the motion of the H/3 and H7 emis- 
sion line peak by measuring their radial velocity. Then 
we constructed periodograms using the AOV algorithm 
(Schwarzenberg-Czerny 1989) and the Scargle (1982) 
method. We scanned a frequency range between 90 c/d 
(the Nyquist frequency) and 1 c/d. Both methods yielded 
consistent results. 

The AOV periodogram is shown in Fig. 4. The main 
peak is found at 16.272 ± 0.058 c/d, along with ± 1 
c/d aliases at 15.3 and 17.3 c/d. Based on the strength 
of the main maximum and assuming a power window 
symmetrical regarding the absolute minimum, we selected 
16.272 c/d (0'?0615 ± 0^0002) as the most likely period, 
corresponding probably to the binary orbital period. Due 
to the possible biased nature of the emission line radial 
velocities, we used the H7 absorption line (see next sec- 
tion) as the best tracer of the white dwarf binary motion, 
measuring the half amplitude of their radial velocity 
which resulted basically independent of the period alias, 
giving confidence about the stellar parameters derived 
in Section 3.5. We found the following ephemeris for the 
inferior conjunction of the secondary star: 



HJD = 2450453.71479(308) + 0.0615(2).B 



(1) 



3.2. Comparison between emission and absorption radial 
velocities 

In this section we compare the radial velocity sets of emis- 
sion and absorption components obtained with different 
methods. It was extremely difficult to find reliable veloci- 
ties for the absorption component. First we tried extract- 
ing the emission core and then fitting a Gaussian func- 
tion to the wings, but the line center was poorly con- 
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strained and we obtained noisy results. Secondly we tried 
deblending the spectral line with emission plus absorp- 
tion Gaussians - using the "d" key in the splot IRAF 
package, obtaining a large amplitude for the absorption 
component (around 110 kms~^). However, numerical sim- 
ulations showed that this result was an artifact: the same 
apparent oscillation was obtained with a synthetic zero 
velocity absorption profile and a moving - 600 km s^^ 
amplitude - emission component (the s-wave, see below). 
Finally, we applied the cross-correlation technique (Tonry 
& Davis 1979), implemented in the fxc IRAF task, to the 
H/3 and H7 white dwarf absorption, excluding the emis- 
sion component inside a range of -2500 km s^^ < v < 
2500 km s~^. Results for H/3 turned to be too noisy, but 
the H7 velocity, binned with the orbital period and fitted 
with a sine function, allowed us to calculate the time of in- 
ferior conjunction of the secondary star. To and the radial 
velocity half-amplitude: K^d = 4 ± 1 kms^^ (Fig. 5). 

In Fig. 5 we observe the contrast between the large 
amplitude of the emission line peak and the white dwarf 
velocity. We will show that the large amplitude of the 
emission peak (619 ± 19 km s~^ for II/3) reflects the char- 
acteristic s-wave observed in the trailed spectra of some 
dwarf novae in quiescence. Note that this s-wave moves in 
antiphase with the white dwarf, consistent with an origin 
close to the inner lagrangian point. 

Fig. 5 also shows that the low K amplitude is incon- 
sistent with noise, as inferred from measurements made 
on synthetic spectra with zero velocity and noise similar 
to the original dataset. We explored the possibility that 
the low velocity pattern in Fig. 5 is the effect of resid- 
ual high velocity emission. However, the s-wave phasing 
and the disk emission phasing (see below) indicate that, 
if present, it should increase the measured amplitude if 
coming from the s-wave and should have minor effects if 
coming from the disk (the effects are summed or canceled 
according to the relative phasing). 

The inferred K^id value is amazing. Their accuracy 
probably lies in the statistical significance of the cross- 
correlation technique, which fits a large sample of pixels si- 
multaneously. Interestingly, low K^d values have also been 
observed in other low mass-transfer rate dwarf nova which 
reveal the white dwarf in their optical spectra, namely WZ 
Sge (40 ± 10 km s'^. Spruit & Rutten 1998), GW Lib 
(40 km s-i, Szkody et al. 2000) and AL Com (0 ± 32 
km s~^, Howell et al. 1998). All the above values were, 
however, obtained from the radial velocity of the emission 
line wing, assumed equal to the white dwarf velocity. On 
the contrary, our results are the first to directly test the 
white dwarf absorption wings with the cross-correlation 
technique. 

For Ha emission, we applied the "double Gaussian" 
convolution mask algorithm (Schneider & Young 1980, 
Shaffer 1983, Home et al. 1986) which provides a robust 
diagnostic method to investigate the behavior of different 
profile sections during the orbital cycle. The method si- 



multaneously shifts two Gaussians of standard deviation 
ag (or alternatively full width at half maximum FWHMg) 
and center separation A along the emission profile until a 
velocity is found for which the convolved flux in both is 
the same. Changing A and FWHMg we can probe differ- 
ent velocity sections of every profile. This method is more 
versatile, allowing us to search for a better trade-off be- 
tween noise and velocity resolution by changing FWHMg. 
It is superior to the simple Gaussian fit for the emission 
component. 

The diagnostic diagram for the Ha emission in Fig. 6 
shows the variations of K and A<I> (the phase shift with 
respect to the ephemeris defined by the absorption com- 
ponent) for velocity sets obtained with different Gaussian 
separations. We find that K is large close to the line center 
but decreases towards the line wing, which could be the 
result of a systematic reduction of the s-wave perturbation 
toward the line wing. This is consistent with the diagnos- 
tic diagram starting with Gaussian separations close to the 
maximum s-wave velocity. In addition, the A<I> behaviour 
indicates that the emission line center precedes the wing 
when crossing the zero phase. Interestingly, the emission 
line wing is almost in phase with the white dwarf mo- 
tion. The above is consistent with wing emission arising 
from the inner disk and following more closely the white 
dwarf motion and low velocity emission coming from a 
hotspot in the disk-stream impact region close to the in- 
ner Lagrangian point. This view is also consistent with 
the fact that K decreases towards the line wing (Fig. 6 
upper panel), maybe approaching the white dwarf value, 
whereas the noise {gk/ K) becomes prominent. 

Phase binned trailed spectra of the continuum and 
white dwarf subtracted spectra (see section 3.5) are shown 
in Fig. 7. These spectra confirm the presence of an addi- 
tional emission component. We find that this component 
is stronger in the higher-order Balmer lines and the He I 
5875 line. This emission component is commonly associ- 
ated with a hotspot in the disk and is called the " s-wave" . 
Disk and hotspot are studied in the next section. 

3.3. Emission line analysis and the optical depth of the 
disk and the hotspot 

We measured the Balmer decrement for both the accre- 
tion disk and the hotspot region to check for different line 
opacities. The Balmer decrement {D{Ha/H(3)) is defined 
as the ratio between the frequency-integrated specific in- 
tensities of the Ha and HP lines. The fact that the disk 
and hotspot emission are simultaneously present in the 
line is a problem for the calculation of decrements. De- 
blending of hotspot and disk emission components is not 
possible without a complex physical model for the disk 
and hotspot regions. We find an approximate solution for 
this problem in the following way. We binned the ob- 
served spectra in 10 phase intervals and corrected them 
for the motion of the white dwarf. The white dwarf spec- 
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trum corresponding to the white dwarf parameters derived 
in section 4 was subtracted from these spectra obtaining 
phase-binned emission profiles. To separate the emission 
components, we assumed symmetrical disk emission and 
separated the total profile for each bin in two components. 
The first consists mostly of disk emission, either the red 
or blue half of the profile depending on the s-wave veloc- 
ity, and a second consisting of disk and hotspot emission. 
The larger and smaller peak had fluxes fmax and fmin re- 
spectively. The hotspot decrement was calculated for each 
phase-binned spectrum assuming that fhotspot = fmax - 
fmin and the disk decrement assuming fdisk = fmin- We 
also assumed a similar shape for the profiles at different 
wavelengths and applied a correction factor to calculate 
the decrements in the fx space. Evidently the above is 
a rough approximation, and phase-bins when the profile 
was almost symmetric (i.e. the hotspot component was 
near the middle of the profile) were not measured. 

Results and comparison with theoretical decrements 
by Williams (1991) are shown in Fig. 8. Williams modeled 
optically thin gas in the emission lines in accretion disks. 
He computed H/3 strengths and Balmcr decrements for a 
grid of temperatures, inclinations and mid-plane accretion 
disk density (A^o in units of nucleons per cm^). Fig. 8 shows 
the results of his models for a representative inclination 
of 52° along with our observations. Results for different 
inclination angles are not significantly different. A direct 
interpretation of Fig. 8 is not possible since the decrements 
are not single- valued functions of disk temperature. How- 
ever, considering some well-established facts, it is possible 
to extract some useful information. First, the observed 
disk decrement {D{Ha/H(3) ~ 2.3) indicates that for any 
realistic disk temperature {Tdisk > few 10'^ K) the Balmer 
lines are optically thin. Second, the hotspot Balmcr decre- 
ment {D(Ha/H/3) « 0.8) is typical of optically thick line 
emission and, if the lines are saturated, it would corre- 
spond to a black-body temperature of 8000iir. The differ- 
ence between hotspot and disk decrement indicates a disk 
optically thinner than the hotspot in the lines. The previ- 
ous analysis was applied to the prototype dwarf nova WZ 
Sge by Mason et al. (2000). These authors found values of 
3.8 and 2 for the disk and hotspot decrements respectively, 
indicating that the disk density in IRXS J105010.3-140431 
is not as low as in WZ Sge. 



ity exponent a (assuming the emissivity ex r~") of an opti- 
cal and geometrically thin Keplerian disk are functions of 
the emission line widths. For Ha we measured 1515, 1915, 
2420 km s~^ at 80, 40 and 10% of the maximum intensity 
level, respectively, yielding R « 0.18 and a w 1.5. The 
emission line widths above were calculated after subtract- 
ing the underlying absorption profile. The corresponding 
figures for the raw emission are 0.22 and 1.6 respectively. 
We observe IRXS J105010.3-140431 close to WZ Sge in 
Fig. 8 by Mennickent & Arenas (1998). The rather large R 
value could indicate inner disk removal. Using equation (4) 
of Mennickent & Arenas (1998) we calculated a supercy- 
cle length of about 5.2 years for IRXS J105010.3-140431, 
which should explain the lack of observed outbursts. This 
recurrence time should correspond to a mass transfer rate 
M - 4 X lO^"* gs"i according to Fig. 2 by Warner (1995). 
This is a very low mass transfer rate similar only to that 
of WZ Sge (M - 1 X 10^^^ gs-i) and HV Vir (M - 3 x 
lO^'' gs-i, Warner 1995). 

3. 5. Synthetic magnitudes and the upper limit for the 
systemic inclination 

Synthetic V magnitudes were calculated from the flux- 
calibrated spectra. Nightly mean values are given in Table 
1. The seeing was stable during the four observing nights 
(less than 0.8 arcseconds), therefore we are confldcnt that 
the slit losses have little effect on the mean magnitudes, 
although they could be the cause of the larger rms com- 
pared with our photometric data. A plot of V versus or- 
bital phase does not reveal eclipses. This confirms the re- 
sults of the differential photometry (Section 2.1). 

We can get an upper limit for the inclination in a 
non-eclipsing CV. For that we used the Egglcton (1983) 
approximation for the volume radius of the Roche lobe of 
the secondary star, obtaining: 



tan(2 imax) 



0A9q 



2/3 



(l-rd/a)[0.6g2/3 + ln(l + (^i/3)] 



(2) 



where q ~ j-p- is the mass ratio and a the binary separa- 
tion. Using the approximation ra/a = 0.6/(1 + g) (Warner 
1995, Eq. 2.61), an upper limit imax ~ 65° is found almost 
independent of q. 



3.4- Evidence for inner disk removal 

In the previous section we have shown that the accretion 
disk of IRXS J105010.3-140431 is similar in some regard 
to that found in long-cycle dwarf novae like WZ Sge. It 
has been suggested that these objects have disks which 
are depleted in their inner regions (e.g. Hameury et al. 
2000, Mennickent 1999, Mennickent & Arenas 1998). One 
method to test this scenario for IRXS J105010.3-140431 is 
looking into the emission line widths. Smak (1981) found 
that the inner to outer disk radius ratio R and the emissiv- 



3.6. Constraints on the stellar masses 

The evolution of cataclysmic variables has been recently 
reviewed with the use of a population synthesis code by 
Howell et al. (2001). The requirement of a Roche lobe- 
filling secondary allows tracing of the evolution of the sec- 
ondary star mass as a function of the orbital period. For 
a given orbital period, two solutions are possible, depend- 
ing on whether the system has P < or P > 0. From 
Fig. 5 of Howell et al. (2001) and using P = 88.6 minutes 
we obtain M2 = 0.14 Mq and M2 = 0.05 Mq respec- 
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tively. The first solution corresponds to a system with a 
nearly main-sequence secondary burning hydrogen in its 
core. The system moves to shorter orbital periods as it 
loses angular momentum by gravitational radiation. The 
second solution corresponds to a system with a secondary 
star having a hydrogen-exhausted degenerated core, like 
a brown dwarf star, and moving to longer orbital periods 
after "bouncing" in the orbital period minimum around 
80 minutes. In order to discriminate between both pos- 
sible scenarios, it is necessary to use a set of empirical 
constraints, which we discuss below. 

The condition that the HWZI of the emission lines 
cannot be greater than the Keplerian velocity correspond- 
ing to the primary radius, combined with the white dwarf 
mass-radius relation by Hamada & Salpeter (1969, yields 
Ml > 0.32 Mq. For this calculation we used the He 5875 
hue width {FWZI = 3680 km s^^), that probably is not 
influenced by underlying absorption. 

The low Kwd also provides strong constraints on the 
stellar masses and systemic inclination through the mass 
function: 



2'kGMi sin^ i = PoKl^ 






(3) 



All the above constraints are plotted in Fig. 9. It is evi- 
dent that a main sequence secondary is incompatible with 
the low K^j^ value. The formal solution gives a secondary 
mass between 10-20 jovian masses for inclination angles 
between 65 and 35 degree. The inclination would have to 
be below 5 degrees for a main sequence companion, which 
is highly unlikely due to the FWZI of He I and the s-wave 
amplitude. Even if for some unknown reason the absorp- 
tion lines do not follow the white dwarf motion, but the 
emission wing velocity tends to this value, the mass func- 
tion for K = 20 km s~^ implies that the inclination should 
be lower than 25 degree for a hydrogen burning compan- 
ion, which is possible but unlikely. The mass transfer rate 



derived in Section 3.4, viz. 6.3 x 10 ^"^Mr^ 



©yr 



and the 



low mass ratio (q < 0.09 derived from the empirical re- 
lationship by Mennickent (1999) using the K^JFWHM 
ratio) nicely fit the prediction of Howell et al. (2001) for a 
post-minimum system with orbital period around 89 min- 
utes. 

We have found the first dynamical evidence for a cata- 
clysmic variable beyond the orbital period minimum. The 
first objects belonging to this class, viz. LL And and EF 
Eri, have recently been discovered by detecting the sec- 
ondary features in the infrared by Howell & Ciardi (2001), 
although indirect evidence has been provided for other ob- 
jects previously, namely AL Com (Howell et al. 1998), WZ 
Sge (Ciardi et al. 1998), V592 Her (van Teeseling et al. 
1999) and EF Eri (Beuermann et al. 2000). 



3. 7. Imaging the accretion disk with Doppler tomography 

The spectroscopic data were phased according to the or- 
bital ephemeris given above. The continuum was sub- 
tracted off in two stages; First, a cubic polynomial was 
fit to the continuum regions between the white dwarf ab- 
sorption dips and subtracted. Second, a 15000 K Pop II 
white dwarf model, convoluted with the instrumental res- 
olution, and with the continuum removed using a poly- 
nomial as above, was fit to the white dwarf absorption 
dips in the data and subtracted to leave only the line 
emission. Using maximum entropy Doppler tomography 
(Marsh & Home 1988) we constructed Doppler maps of 
the Ha, H/3, H7, US and He I 5876 lines. The data were 
not phase-binned before constructing the Doppler maps; 
this maximized phase coverage with 57 spectra covering 
all orbital phases observed. The maps are shown in Fig. 10. 
The Ha map shows a symmetric disk and a weak hotspot 
with velocities between that of the ballistic stream and 
local Keplerian velocity. In the H/?, H7, US maps the disk 
becomes progressively fainter, while the hotspot becomes 
stronger and shifts to velocities closer to the local Kep- 
lerian velocity. This is consistent with high order Balmer 
lines formed deeper in the ballistic stream, where the ex- 
cess of kinetic energy is radiated by interaction with the 
disk. The He I 5876 map shows low velocity bright spot 
emission; this may indicate substantial photoionization of 
the gas stream and/or that the He I bright spot originates 
further up the ballistic stream than the HI bright spot, 
suggesting that the accretion disk extends further out than 
would be deduced from the H I emission alone. In the nova- 
like system V347 Pup for instance, the Hell emission is 
produced at inner parts of the disk when compared with 
HI (Diaz & Hubeny 1999). 

The Ratiocd Doppler Map for IRXS J105010.3-140431 
is shown in the bottom right of Fig. 10; this was con- 
structed by dividing the Ha map by the H/3 map shown 
in Fig. 10. The change in the Balmer decrement in the 
region of the bright spot indicates hotter optically-thick 
emission than is seen from the disk. It is interesting to note 
the similarity in the structure of the ratioed Doppler maps 
of 1RXSJ105010.3-140431 and that of WZ Sge (Skid- 
more et al. 2000). However, the ratios of the disk and 
hot spot emission measured in IRXS J105010. 3-140431 
are generally lower than measured in WZ Sge, indicat- 
ing higher temperatures and/or greater optical depth in 
IRXS J105010.3-140431, as found in Sec. 3.3. 



3.8. Radial brightness temperature profile 

Using the same method as described in Skidmore et al. 
(2000), the radial profile of the ratioed Doppler map was 
calculated by finding the average ratio at each radius in 
velocity space. Areas of the ratio map affected by the 
bright spot flux were masked when calculating the radial 
profile. These areas covered between -20 and 75 degrees 



Mennickent et al.: Discovery of a CV with a sub-stellar companion 



measured anti-clockwise from the y axis. The mean line 
ratio at each velocity was used to determine a tempera- 
ture assuming a simple blackbody source function at the 
wavelength of the line emission. This oversimplification 
is used as a rough guide to infer temperature gradients 
in the disk. The system parameters for IRXS J105010.3- 
140431 determined in this paper were used to convert the 
velocity space radial temperature profile into a real space 
radial temperature profile. The radial temperature profile 
for IRXS J105010.3-140431 is shown in Fig. 11. 

We note that the inner disk temperatures are larger 
than those suggested by the simple decrement analysis. 
Apparently the decrement study gives more weight to the 
outer disk regions. This is easily explained since there is 
a larger contribution of these low velocity regions to the 
emission profile and the decrement study is based on the 
integrated flux along the profile. 

As the slope of the radial temperature in the disk is 
close to T ^ r""^/**, we can see that the disk is in a steady 
state. In quiescent dwarf novae the disk radial temperature 
profile is usually observed to depart from the r^'^l^ law 
and to have a much flatter profile (Wood, Home & Vennes 
1992, Skidmore et al. 2000). The departure from T ~ r'^/^ 
occurs at radii of r^n ^0.1 and ro„j ^0.5. If these radii 
represent the inner and outer edges of the disk then R= 
^^^^ ~0.2 as found in Sec. 3.4. 



4. Modeling the absorption spectrum 

The observed spectrum of IRXS J1G5G1G.3-140431 clearly 
shows broad Balmer absorption profiles flanking the emis- 
sion lines. These features present extended line wings typ- 
ical of Stark broadened Balmer lines; therefore, we may 
associate them with the underlying white dwarf spectrum, 
instead of Doppler-broadened absorption from the accre- 
tion disk photosphere. The conspicuous presence of the 
Balmer absorption component indicates that the white 
dwarf photosphere is bright when compared to the accre- 
tion disk continuum. In most accreting cataclysmic vari- 
able systems the accretion disk flux prevents the obser- 
vation of the compact star. The detection of the white 
dwarf in IRXS J105010.3-140431 offers the opportunity 
of obtaining basic photospheric parameters for the white 
dwarf by comparing the observed profiles with a grid 
of synthetic spectra. On the other hand, helium absorp- 
tion features could not be identified in the profiles of 
He IAA6678, 5876, 4471 A. This fact suggests a DA clas- 
sification for the white dwarf. 



J^.l. High-gravity atmosphere models and spectrum 
synthesis 

In order to obtain a grid of model DA white dwarf spectra, 
pure hydrogen atmospheres were calculated assuming hy- 
drostatic and radiative equilibrium. Only the continuum 
and hydrogen opacity are included, along with scattering 



and free particle processes. The plane-parallel atmosphere 
structure is solved locally at 70 mass-depth points by the 
complete linearization method using the code TLUSTY 
(Hubeny & Lanz 1995). 

Once the density, temperature and opacity is calcu- 
lated, the radiative transfer equation is solved for each 
frequency in the region of interest by the program SYN- 
SPEC (Hubeny, Lanz & Jefferys 1994). NLTE departures 
are allowed for the hydrogen levels population. The im- 
proved calculations of the Stark broadening by Shoning & 
Butler (1989) in addition to the thermal Doppler compo- 
nent arc used to compute the Balmer line profiles. Both 
emergent flux {H\) and theoretical continuum are com- 
puted in absolute units. Our line profile calculations were 
found to be in good agreement with the white dwarf mod- 
els by Jordan & Kocster (1986). 

4-.2. Model fitting and white dwarf parameters 

About 90 model spectra were computed with effective 
temperatures ranging from 10000 K to 65000 K and \ogg 
between 6.0 and 9.5. Continuum subtracted and contin- 
uum normalized model sets were prepared using the the- 
oretical continuum calculated in the absence of bound- 
bound transitions. Subsequently we fitted a low order 
polynomial pseudo-continuum to the models in the region 
covered by our spectroscopy. This residual continuum cor- 
rection of model spectra is necessary for their comparison 
with the observations because of the significant overlap- 
ping of the Balmer line wings beyond H/3. 

The regions for fitting were carefully selected in the 
observed spectrum to exclude the contribution from the 
narrow central emission and from weak emission lines in 
the extended wings. Two different procedures were em- 
ployed to fit model spectra to the observations: 

The first method consists of fitting the line profile 
in flux units using continuum-subtracted spectra. Each 
model spectrum in the grid is arbitrarily scaled to match 
the integrated flux in the observed absorption lines and 
then the reduced chi-square statistic is computed for each 
pair (Teff, log 5 ). This method does not require a pre- 
vious knowledge of accretion disk continuum contribution 
to the total fiux. However, it has the disadvantage of dis- 
regarding the line to continuum fiux ratio. Results from 
this method are shown in Fig. 12. A valley in the x^ sur- 
face is evident from the graph, containing two minima. A 
sample individual fit can be seen in Fig. 13. 

A caveat related to the effect of the unknown interstel- 
lar extinction should be made here. This effect, in prin- 
ciple, must be considered when the simultaneous fitting 
of Ha, H/3 and H7 is performed. However, a comparison 
with the fit of H/3 and H7 alone does not show significant 
discrepancies. In addition, the expected extinction at the 
target galactic latitude (6/7 = 39 deg) should be low. 

The second method uses the continuum-normalized 
spectrum as input and deals with the veiling by the other 
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continuum sources in the binary. Unfortunately, the shape 
of the underlying disk continuum emission is unknown. 
However, given its low luminosity relative to the white 
dwarf one may expect under standard theory that such a 
disk should present relatively low temperatures. For the 
sake of simplicity it was assumed that the disk continuum 
shape is similar to the white dwarf continuum, but scaled 
by an unknown factor. The relative disk contribution to 
the total flux (/) is independent of the wavelength and 
may be left as a free parameter in the model fitting. Typ- 
ical results are shown in Fig. 14, where the behavior is 
similar to that seen in Fig. 12. Simulations show that the 
minima drift toward lower temperatures and higher grav- 
ities as the veiling factor / increases. The increase in / 
also implies a significant decrease in the quality of the fit. 
It became clear from our simulations that the disk should 
contribute less than 20% of the system flux. 

The dynamical solution discussed in Section 3.6 gives 
a range for the white dwarf mass too large to discriminate 
between the two minima found in the x^ plots. By com- 
bining the quality criteria obtained in both procedures we 
estimate log g ranging from 7.6 to 9.3 and T^ff from 13 000 
K to 24 000 K. 



5. Conclusions 

— A spectroscopic analysis of the ROSAT target 
1RXSJ105010.3-140431 has been realized noting the 
close similarity between this source and the prototype 
dwarf nova WZ Sge. In particular, the white dwarf 
photospheric lines are revealed in the spectrum, which 
is a signature of a very low mass accretion rate; 

— We find a most likely orbital period of 88.6 minutes 
although ± 1 cycle/day aliases cannot be ruled out; 

— The absorption spectrum can be modeled with a DA 
white dwarf with 13 000 K < T^ff < 24 000 K; 

— The mass function suggests that the system is beyond 
the orbital period minimum having an undermassive 
secondary; 

— As in other short orbital period dwarf novae systems, 
evidence for inner disk removal is found. In addition, 
using empirical relationships for dwarf novae, we pre- 
dict a supercycle of several years for this object; 

— The Doppler maps and the Balmer decrement study re- 
veal significant differences between the optical depth 
of the disk and hotspot. We find an optically thick 
hotspot with T ~ 8000 K and an optically thinner 
disk. The hotspot has velocities consistent with dissi- 
pation of kinetic energy in the stream-disk interacting 
region, and it is the main emitting source in high order 
Balmer lines and the He I 5875, but not in Ha, where 
the emission is mainly distributed along a symmetri- 
cal disk. Contrary to the Balmer lines, the helium line 
seems to be formed in the outer stream-disk impact 
region; 



— We constructed the temperature profile for the disk 
and, contrary to that found in other dwarf novae, we 
find an almost steady state standard disk, with T ~ 
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Figure Captions 

— Fig.l An ESO Digitised Sky Survey image indicating 
the comparison stars used for differential photometry. 
Tlie field of view is 5'x 5'. 

— Fig. 2 The differential magnitudes in the V band along 
with their error bars. 

— Fig. 3 The averaged spectrum (upper panel) and a close 
up view of the main spectral lines (below panel). From 
up to down the He 5875, H7, H/3 and Ha lines are 
shown. 

— Fig. 4 The AOV periodogram of the H/3 emission line. 
The maximum at 16.27 c/d ((y?0615) is identified as 
the orbital period. 

— Fig. 5 Upper panel: The H/3 (squares) and H7 (circles) 
s-wave velocity folded with the ephemeris given by 
Eq. 1. Below panel: The H7 absorption line velocity 
(squares) obtained by cross correlation and the best 
sine fit along with results for zero velocity synthetic 
absorption profiles (circles). 

— Fig. 6 The diagnostic diagram for the Ha emission line 
considering the ephemeris given by Eq. (1). 

— Fig. 7 The trailed spectra of the emission lines. The 
data have been binned into 20 phase bins and are re- 
peated for clarity. Note how the relative strength of the 
bright spot emission with respect to the disk emission 
increases from Ha to US. A double peak is evident in 
Ha. 

— Fig. 8 Theoretical Balmer decrements for optically thin 
disks emission lines for inclination of 52° according to 
Williams (1991). The models for T= 8000 K, 10 000 
K, and 15 000 K, all approach to D(Ha/H/3) = 1 at 
high densities, as expected for optically thick lines. The 
lines indicate orbit-averaged decrements for disk and 
hotspot and their corresponding uncertainties. 

— Fig. 9 The Mi — M2 plane showing the dynamic solu- 
tions for Ki = 4 km s^^ (absorption wings) and Ki 
= 20 km s~^ (emission wings). The vertical lines in- 
dicate possible secondary star masses for the orbital 
period accordingly to results of the evolutionary code 
by Howell et al. (2001). A post orbital period mini- 
mum secondary is favoured by our velocities. 

— Fig. 10 Doppler maps of the Ha, H/3, H7, US and He I 
5876 emission lines. Overlayed is the Roche geometry 
of the system with Mi = 1.0 Msun,M2 = 0.05Ms„„, 



Porb — 0.0615 d and i — 65 deg. The dashed line shows 
the local Keplerian velocity along the stream path as- 
suming circular orbits around the primary. The cir- 
cular line represents 1/3 of the primary Roche lobe 
radius. The relative increase in bright spot strength 
with respect to the disk emission is evident as we go 
from Ha to US. Note that the He I bright spot is posi- 
tioned with a slower velocity than the H I bright spots. 
In general the bright spots have velocity between that 
of ballistic stream and the local Kepler velocity. Bot- 
tom Right: The Ratio map of the Ha and H/3 Doppler 
maps. Note the striking similarity in the structure of 
the ratio map of IRXS J105010.3-140431 to that of WZ 
Sge shown in Fig. 6 of Skidmore et al. (2000). 
Fig. 11 The radial temperature profile for the disk emis- 
sion in IRXS J105010.3-140431 assuming a blackbody 
source function. The diagonal dashed line indicates the 
slope of the radial temperature profile for a standard 
alpha disk where T ^ r~^'^. 

Fig. 12 x^ as a function of surface gravity and temper- 
ature according to the results of Method 1 . 
Fig. 13 Illustration of the spectrum synthesis analysis. 
Fig. 14 x^ as a function of surface gravity and temper- 
ature according to the results of method 2. The plot 
corresponds to a value of / = 0.08 (see text). 
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